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FOREWORD 


This  report  presents  the  results  of  parametric  studies  on  several 
families  of  rocket  grain  configurations  obtained  by  photoelastic 
analysis.  These  results  are  presented  graphically  in  the  form  of  a 
stress  factor  H  as  a  function  of  the  various  parameters.  In  addition, 
empirical  relationships  have  been  derived  for  most  of  the  studies 
conducted  which  represent  the  observed  dependency.  A  section  on  the 
application  of  the  test  results  to  more  general  engineering  problems 
is  presented.  In  this  section  a  comparison  is  made  between  a 
recommended  method  of  application  and  a  limited  number  of  numerical 
solutions  obtained  via  computer  analysis. 

The  authors  wish  to  acknowledge  the  able  assistance  in  the  conduct 
of  the  experimental  portion  of  this  work  of  A.  P.  Waggoner  and  W.  L. 
Fourney . 
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ABSTRACT 


This  report  presents  the  results  of  parametric  studies  on  several 
families  of  rocket  grain  configurations  obtained  by  photoelastic  analysis. 
In  an  earlier  paper  the  authors  described  ?  parametric  investigation 
of  a  family  of  grains  which  were  defined  by  three  parameters:  N  the 
number  of  star  points,  a/b  the  port  fraction,  and  a/p  the  fillet 
radius  factor.  The  results  of  that  study  led  to  the  following  empirical 
formula  for  the  stress  concentration  factor, 

»-i/3  [i ♦*!/*} 

In  the  present  investigation  the  study  has  been  extended  to  include 

families  requiring  four  parameters  to  complete  their  description.  Four 

additional  families  have  been  investigated  tc  determine  the  effect  of 

a)  slot  width,  b)  positive  wedge  angles  of  the  star  slot,  c)  negative 

wedge  angles  of  the  star  slot,  and  d)  eccentricity  of  elliptically 

-1/3 

shaped  star  tips.  In  each  of  the  studies  It  is  shown  that  the  N  ' 
rule  holds  approximately,  and  empirical  formulas  similar  to  the  one 
given  above  have  been  derived  in  some  cases 

A  section  on  the  application  of  the  test  results  to  more  general 
engineering  problems  is  presented.  In  this  section  a  comparison  is 
made  between  a  recommended  method  of  application  and  a  limited  number 
of  numerical  solutions  obtained  via  computer  analysis.  Under  the 
worst  set  of  conditions  available,  the  streus  calculations  agreed 
within  101  and  the  strain  within  15%. 

As  a  result  of  the  large  number  of  parameters  used  in  this  study, 
it  is  practical  to  obtain  the  maximum  stress  for  most  grain  configura¬ 
tions  using  the  data  presented. 
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I. 

INTRODUCTION 

In  the  design  of  solid  propellant  rocket  engines,  the  use  of 
internal  slot  configurations  has  complicated  the  analysis  of  stress 
and  strain.  Additional  complications  enter  into  the  problem  because 
of  the  use  of  viscoelastic  materials,  elastic  cases,  and  tapered 
shapes  and  because  of  end  effects  and  finite  deformations.  Ignoring 
these  complications  for  a  moment,  we  might  consider  the  problem  of  an 
infinitely  long,  star-perforated  cylindrical  body  of  linearly  elastic 
material  subjected  to  infinitesimal  deformations.  Even  the  solution 
of  this  highly  simplified  problem  defies  exact  analytical  an? lysis  1 
However,  it  may  be  solved  by  approximate  methods  of  analysis  or  by 
using  the  techniques  of  photoelasticity .  If  the  solution  to  this 
problem  has  some  application  in  the  solution  of  the  even  more  diffi¬ 
cult  finite  len^ili,  case-bonded,  viscoelastic  problem,  then  it  serves 
a  worthwhile  purpose.  It  will  be  shown  that  this  is  indeed  the  case. 

In  view  of  this,  several  photoelastic  investigations  of  stresses 

(1  2  3)* 

in  solid  propellant  rocket  grains  have  been  conducted,  ’  ’  This 
report  presents  the  results  of  a  series  of  parametric  tests  on  four 
families  of  grain  configurations  The  geometry  of  each  family  is 
completely  characterized  by  four  parameters. 

The  present  results  are  compared  to  previous  work  by  the 
authors  the  agreement  is  excellent.  An  error  analysis  of  the 
data  has  been  conducted;  the  method  of  analysis  and  the  results  are 
outlined.  The  analysis  indicates  that  the  data  which  is  presented  is 
accurate  within  +  5%  over  the  range  of  interest. 

^Numbers  enclosed  in  parentheses  refer  to  References  listed  on  Pa^e  48. 
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II. 


GENERAL  THEORY 


The  solution  to  the  problem  of  en  infinitely  long,  star-perforated 
cylindrical  body  (Figure  1-a)  is  (by  symmetry)  independent  of  the 
'xial  coordinate  (which  is  taken  to  be  the  z  coordinate),  provided 
the  loading  on  the  yli.idri  cal  surfaces  is  independent  of  z  and  the 
z  corap.-ent  of  all  loads  is  zero.  The  problem  is  classified  as  a 
problem  in  plane  -train.  If  the  material  is  homogeneous,  isotropic, 
and  linearly  elastic  and  if  the  body  forces  are  zero,  then  the  follow¬ 
ing  equations  can  be  shown  to  govern  the  problem. ^ 
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Constitutive  equations  relating  strain  to  stress 
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Auxiliary  equations  for  calculating  the  remaining  stress 
and  strain  components 
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Now  consider  a  thin  slice  cut  from  this  cylinder  and  subjected  to  the 

same  load  per  unit  area  on  the  cylindrical  surfaces  (see  Figure  1-b). 

The  upper  and  lower  surfaces  are  assumed  to  be  free  of  tractions. 

This  problem  is  classified  as  a  problem  in  plane  stress,  and  its 

(4) 

solution  is  governed  by  the  following  set  of  equations. 


Field  equations  for  the  stress  components  Oxx,  CTyy,  (7^ 
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Constitutive  equations  relating  stress  to  strain 
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Auxiliary  equations  for  calculating  the  remaining  stress 
and  strain  components 


(6) 


In  each  problem,  the  field  equations  are  to  be  solved  subject  to 
certain  boundary  conditions.  The  stress  components  0^x,  ^yy*  ^xy 
obtained  from  the  solution  are  then  substituted  into  the  remaining 
equations  to  find  the  complete  stress-strain  field.  Notice  that  the 
field  equations  of  plane  stress  and  plane  strain,  with  zero  body  force, 
are  identical.  Thus,  if  the  stress  boundary  conditions  for  a  problem 
in  plane  strain  are  identical  to  the  stress  boundary  conditions  for  a 
problem  in  plane  stress,  the  two  solutions  for  ^xx.  °yy»  (“Jxy  must 
be  the  same  and  will  be  independent  of  material  properties  as  long  as 
E  and  V  do  not  enter  the  boundary  conditions.  This  observation  is 
the  basis  for  the  application  of  photoelastic  analysis  to  the  rocket 
grain  (plane  strain)  problem. 

Let  us  consider  the  specific  case  of  a  long  circular  cylinder 
with  a  star-shaped  cylindrical  perforation  (Figure  1-a).  If  there 


are  no  body  forces  and  if  Che  surface  tractions  are  a  uniformly  dis¬ 
tributed  pressure  pQ  on  the  external  cylindrical  surface  and  a  uni¬ 
formly  distributed  pressure  p^  on  the  internal  cylindrical  surface, 
then  the  boundary  conditions  are: 


External  Surface: 

Outward  Normal  n 

— ► 

Surface  Traction  T  “  -p  n 

o 

Boundary  Conditions 


°ixnx  +  C^xyny  -  -Po^ 
°xynx  +  ^yyVy  m  'Po^ 


Internal  Surface: 

Outward  Normal  m 
Surface  Traction  T  =  -pjm 
Boundary  Conditions 

°"x*“x  +  ^xy^y  ”  "Pi11^ 

V*  +  (^yyIDy  "  -Pi°y 
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(8) 


These  boundary  conditions  are  independent  of  E  and  V  .  Therefore, 
a  geometrically  similar  thin  disk  of  any  linearly  elastic  material, 
loaded  with  the  same  internal  and  external  pressure,  will  be  subjected 
to  the  same  stress  components  CT^,  0~yy,  .  Hence  the  stress 


field  in  the 


The  photoelastic  test  is  performed  in  the  following  manner.  A 

thin  disk  model  of  the  appropriate  shape  is  prepared  from  a  material 

that  exhibits  birefringence  when  subjected  to  strain.  This  model 

could  then  be  loaded  with  internal  pressure  and  external  pressure 

p'  and  viewed  in  a  circular  polariscope.  The  isochromatic  fringe 
o 

pattern  which  would  be  observed  is  directly  related  to  the  difference 
in  principal  stress,  -  (J^  ,  at  each  point  in  the  model. 

In  practice,  it  is  difficult  to  apply  a  uniform  pressure  p^  to 
the  geometrically  complicated  inner  boundary.  It  is  much  more  con¬ 
venient  to  apply  only  p^  to  the  simple  circular  outer  boundary. 
Fortunately  the  solution  for  internal  and  external  pressure  can  be 
easily  deduced  from  observations  made  on  a  model  with  external  pressure 
only.  This  result  follows  from  the  well  known  principle  of  super¬ 
position  for  the  linear  theory  of  elasticity. 

Consider  two  loading  conditions  applied  to  the  disk.  Loading  A 
consists  of  equal  internal  and  external  pressure  pj^  .  Loading  B 
consists  of  external  pressure  p^  -  p[  and  zero  internal  pressure. 

It  is  clear  that  the  superposition  of  the  solutions  to  theLe  two 
problems  is  the  solution  for  the  problem  with  internal  pressure  p[ 
and  external  pressure  p^  .  Problem  B  is  convenient  for  photo¬ 
elastic  analysis.  Problem  A  is  a  case  of  two-dimensional  hydrostatic 
loading  and  has  the  solution 


(9) 


All  directions  in  the  x,y  plane  are  principal  directions.  Thus 

-  CT2  -  -p'  (9-a) 

Since  the  equations  which  govern  problem  B  are  linear,  it  is 
clear  that  the  magnitude  of  any  stress  component  at  any  point  in  the 
disk  is  linearly  proportional  to  fie  magnitude  of  the  external  pressure. 
If  we  let  p0  *  Pq  -  p|  and  if  we  designate  the  magnitude  of  the  maxi¬ 
mum  stress  which  occurs  in  the  solution  to  B  for  a  given  value  of  pQ 
by  crM  ,  then  CTj^/pQ  will  be  a  constant  for  all  values  of  pQ 
The  determination  of  the  value  of  this  constant,  as  a  function  of 
various  parameters  describing  the  star  geometry,  is  the  primary  objec¬ 
tive  of  this  study. 

In  star  configurations  loaded  with  external  pressure  r>o  ,  the 
maximum  stress  is  found  to  occur  on  the  inner  boundary  at  the  star 
tip.  Since  the  inner  boundary  is  free  from  tractions,  the  principal 
directions  along  the  boundary  are  the  directions  perpendicular  to  and 
tangent  to  the  boundary.  The  principal  stress  acting  on  the  surface 
with  a  normal  perpendicular  to  the  boundary  is  aero.  The  principal 
stresr  acting  on  a  surface  with  a  normal  parallel  to  the  boundary,  at 
the  point  on  the  boundary  where  the  maximum  fringe  order  (difference 
in  principal  stress)  occurs  is  (7^  .  Since  one  principal  stress  is 

zero,  it  follows  that  the  magnitude  of  the  difference  of  principal 
atresa  j -  Oj  J  "  |  0#  J  .  The  difference  in  principal  streaa  la 
directly  related  to  the  fringe  order  n  observed  at  the  point  by  the 


(10) 
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where  K  Is  the  fringe  constant  and  t  Is  the  model  thickness;  hence, 
can  be  directly  determined  from  observation  of  the  fringe  order  n 
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III. 


TESTING  PROCEDURES 

A  test  jig  has  been  designed  which  allows  a  uniform  pressure  pQ 
to  be  applied  to  the  external  periphery  of  the  model.  In  practice 
this  pressure  is  supplied  by  a  regulator  from  a  high  pressure  nitrogen 
bottle  and  may  be  varied  between  0  and  500  psi.  The  ratio 
can  be  obtained  from  observations  at  one  pressure.  However,  random 
error  will  be  reduced  if  a  series  of  observations  are  made  at  various 
press  ires.  The  resulting  curve  of  CT^  versus  pq  should  be  a  straight 
line,  end  its  slope  is  the  desired  ratio  (XM/p 

The  direct  observable  in  these  experiments  is  the  fringe  order  n  , 
which  has  been  shown  to  be  proportional  to  CT^  .  The  experimental 
procedure  therefore  is  to  measure  the  pressure  p^  which  is  required 
to  cause  a  first  order  fringe  at  the  concentration  point,  the  pressure 
p^  required  to  cause  a  second  order  fringe,  etc.,  up  to  some  maximum 
fringe  order  dictated  by  the  pressure  limitations  of  the  system  or  by 
the  limits  of  elasticity  of  the  model.  For  the  present  series  of  tests 
the  pressure  limit  was  usually  reached  between  the  sixth  and  eighth 
fringe  order.  Thus  six  to  eight  data  points  were  available  on  the  p^ 
versus  n  curve,  and  a  precise  slope  could  be  easily  determined. 
Moreover,  either  four  or  eight  equal  stress  concentrations  appear  in 
the  model  (The  number  depends  upon  geometry.  See  Figures  2  and  3.) 
so  that  this  procedure  can  be  repeated  four  or  eight  times,  and  a 
weighted  average  of  the  slopes  is  then  used  to  determine  the  best  value 
of  <XM/p0  for  the  particular  geometry  being  tested.  A  least  squares 
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data  reduction  method  has  been  used.  It  is  described  later  in  the 
report.  A  typical  set  of  curves  of  pressure  versus  fringe  order  is 
shown  in  Figure  4. 

In  order  to  convert  the  direct  experimental  result  n/pQ  (fringe 
order /’pressure)  inL:>  the  desired  quantity  (7m/ Po  »  is  necessary 
to  know  the  fringe  constant  K  in  the  stress  optic  law,  equation  (10). 

The  fringe  constant  has  been  determined  for  each  sheet  of  material 
(CR-39)  used  in  these  tests.  The  value  was  derived  from  several  tests 
in  which  the  fringe  order  versus  the  load  F  was  observed  in  a  specimen 
for  which  (Tj  -  02  versus  load  was  known  from  analytical  considera¬ 
tions.  The  particular  test  which  was  used  was  the  diametral  compression 

(  a) 

of  a  solid  disk  (Figure  5).  The  stress  solution  of  this  problem  isx 
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A  typical  curve  of  fringe  order  versus  F  is  shown  in  Figure  5.  As  an 
example  of  the  consistency  of  the  tests,  ten  specimens  were  tested  from 
Sheet  #1,  and  the  range  of  K  observed  was  93.4  +  1.0  psi-in/fringe. 

Two  additional  tests  were  used  to  check  the  overall  validity  of 
the  test  procedure.  In  one  test,  10  thick-walled  cylinders  were  loaded 
with  external  pressure  in  the  pressure  jig.  Since  the  stress  solution 
is  known  in  this  case  from  theory,  the  correlation  between  test  results 
and  theory  can  be  observed.  The  test  results  are  shown  in  Figure  6. 
Unfortunately,  the  test  is  not  too  sensitive  since  the  stress  levels 
for  a  given  amount  of  external  pressure  are  low,  so  that  within  the 
500  psi  test  limits  only  3  or  4  fringes  can  be  observed.  In  addition, 
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the  stress  gradients  are  low  so  that  the  fringes  are  not  sharply  defined. 
Nevertheless ,  the  agreement  between  theory  and  experiment  is  good  and, 
because  of  the  limitations  inherent  in  the  test,  this  discrepancy 
probably  indicates  the  maximum  overall  deviation  to  be  expected  between 
the  true  values  of  H  and  the  experimental  values  in  this  type  of 
testing. 

To  check  consistency  with  previous  experimental  work,  a  number  of 
simple  slot  configurations  were  tested  (Figure  7-a).  The  result  of  the 
present  tests  is  compared  to  previous  test  results  in  Figure  8.  Para¬ 
metric  studies  of  this  configuration  have  been  published  in  Reference  3. 
Excellent  agreement  was  found  between  the  two  sets  of  data. 


IV. 


DATA  REDUCTION 


Typical  data  from  a  test  on  one  specific  geometry  consists  of 
the  pressures  corresponding  to  the  first  6  to  8  fringe  orders  for  each 
of  the  locations  where  the  maximum  stress  occurs.  The  objective  of 
the  data  reduction  process  is  to  extract  from  this  set  of  50  or  60 
values  of  pressure  the  best  value  for  the  ratio  -CX^/po  • 

As  a  first  step,  the  best  value  of  ^m/pq  *s  determined  for 
each  of  the  locations  where  maximum  stress  occurs.  0^/po  Is  related 
via  the  fringe  constant  to  the  slope  of  the  best  straight  line  through 
the  data  presented  in  a  pressure  versus  fringe  order  plot.  If  all  of 
the  data  is  given  equal  weight,  the  slope  of  the  best  straight  line 

(5)  T  £ 

through  the  data  can  be  found  by  a  least  mean  squares  process.  If 
the  data  is  thought  of  as  a  set  of  M  pairs  of  numbers  (x,*y,)  , 

then  the  best  slope  can  be  shown  to  be 


A 


M 

M^xyri 


M 

Mix, 

i-l 


2 


M 

(I*,) 

t-1 


2 


(12) 


where 


y  *  AX  +  B  .  The  standard  deviation  of  A  can  be  shown 


to  be 


03) 


12 


In  this  way,  L  values  of  best  slope  A£  ,  with  corresponding  stand¬ 
ard  deviations  Sj  ,  are  calculated.  The  best  value  of  A  is  then 
taken  to  be  a  weighted  average  of  these  L  values. 


A  ■  best  estimate  of  A 


ki,s* 

X  1/s/ 

j«i  J 


The  standard  deviation  of  A  is  calculated  from 


(14)  j 

* 

j 
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S_ 

A 


(L  -  1)  Xl/S, 

J-l 


(15) 


DATA  REDUCTION  EXAMPLE 


The  following  raw  data  was  obtained  from  Test  #29 


Table  I 


Fringe 

Order 


External 
s.  Pressure 


Location  of  Stress  Concentration 

#1 

#2 

#3 

#4 

#5 

#6 

#7 

#8 

62 

58 

61 

55 

64 

62 

66 

65 

115 

105 

113 

111 

143 

113 

115 

122 

179 

169 

171 

172 

177 

173 

179 

187 

231 

223 

227 

?34 

228 

229 

230 

244 

293 

284 

284 

295 

293 

284 

288 

305 

357 

347 

348 

360 

357 

350 

352 

373 

6 .44 

6.52 

6.63 

6.21 

6.45 

6.61 

6.62 

6.18 

0.10 

0.12 

0.09 

0.07 

0.14 

0.10 

0.12 

0.07 

The  best  slope  A  is  calculated  for  each  point  of  stress  concen¬ 
tration  using  equation  (12).  For  exaiuple,  fc~  stress  concentration  #1, 
M  •  6,  t  »  0.246,  A  ■  6.44,  and  S  •  0.10.  The  value  of  A  for  each  of 
the  eight  points  of  stress  concentration  is  shown  in  the  above  table, 
with  corresponding  values  of  S  equation  (14)  is  now  used  to  find 
A  the  best  value  of  A  .  with  L  ■  8 


Equation  (15)  gives  the  standard  deviation  of  A  . 


S_  =  0.07 
A 

Thus,  for  the  geometry  of  Test  #29,  the  bracket  H  ■  6.40  +  0.07  con¬ 
tains  the  correct  value  of  H  with  a  probability  of  677o.  The  bracket 
H  ■  6.40  +  0.14  contains  the  correct  value  of  H  with  a  probability 
of  95%. 
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VI. 


TEST  RESULTS 


The  tests  conducted  during  this  program  have  established  the 

stress  variations  which  arise  when  the  simple  slot  configuration  shown 

in  Figure  7-a  is  modified ,by  either  changing  the  shape  of  the  tip  or 

the  slope  of  the  slot  wails.  An  extensive  series  of  tests  of  the 

(3) 

simple  slot  configuration  has  been  previously  reported.  For  com¬ 
pleteness,  the  results  are  shown  in  Figures  9  through  14.  The  simple 
slot  configuration  is  characterized  by  three  parameters:  a/b  ,  a/p  , 
and  N  .  The  variation  of  H  as  a  function  of  a/p  ,  with  N  and 
a/b  fixed  is  shown  in  Figures  9  through  13.  A  typical  photoelastic 
fringe  pattern  is  shown  in  Figure  2.  One  of  the  more  interesting 

results  of  this  study  was  that  the  variation  of  H  as  a  function  of 

-1/3 

N  ,  with  a/p  and  a/b  fixed,  was  of  the  form  H  ~  N  '  over 

the  range  tested,  except  near  the  limit  points  a/p  -*»  1  .  (It  is 
clear  that  for  all  values  of  N  and  p  ,  when  a  ■  p  ,  the  simple 
8 lot  configuration  becomes  a  circular  port.)  This  behavior  is  shown 
in  crossplots  (Figure  14).  The  following  empirical  formula  was  derived 
for  H 


(16) 


and  as  shown  in  Figures  9  through  14  fits  the  data  within  5%  for 
a/p  >  4  . 
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In  the  present  series  of  tests  the  simple  slot  configuration  was 
modified  by 

a)  Widening  the  slot  tips — see  Figure  15. 

b)  Replacing  the  parallel  slot  walls  with  nonparallel 
walls--both  positive  and  negative  angle--see  Figures 
16  and  17. 

c)  Replacing  the  semicircular  star  tip  with  a  semiellipse — 
see  Figure  18. 

The  following  sections  of  the  report  give  further  details  on  the 
configurations  tested  and  the  values  of  H  observed. 
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VII. 


SLOT  WIDTH  FACTOR  TESTS 

A  series  of  tests  was  run  to  determine  the  effect  of  slot  width 
on  maximum  stress.  The  geometry  is  shown  in  Figure  15.  It  is  completely 
defined  by  the  four  parameters:  N  the  number  of  star  points,  a/b 
the  port  fraction,  a/p  fillet  radius  factor,  and  d/2p  the  slot 
width  factor.  The  majority  of  the  tests  were  of  configurations  with 
N  ”  4  .  In  addition,  6  tests  were  made  of  configurations  with 
N  f  4  .  The  results  of  the  tests  with  N  *  4  are  presented  in 
Figures  19  through  22  by  holding  d/2 p  constant  at  integral  values 
from  2  to  5  inclusively  and  plotting  H  versus  a/p  for  constant 
values  of  a/b  . 

The  slot  width  factor  is  varied  between  the  limit  point,  where 
d/2 p  ■  1  ,  which  represents  a  simple  slotted  grain  (See  Figure  8), 

and  the  limit  point  where  d  is  allowed  to  increase  to  the  value 
d  ■  2a  .  These  limiting  geometries  are  shown  in  Figure  7.  Both  of 
these  limit  points  were  established  by  means  of  photoelastic  tests. 

The  limit  point  corresponding  to  d/2 p  *1  is  identical  to  the  con¬ 
figuration  treated  in  Reference  3;  however,  17  additional  tests  were 
conducted  to  confirm  the  correlation  between  these  two  series  of  tests. 
The  other  limit  point  was  established  by  photoelastic  tests  on  the 
geometry  shown  in  Figure  7-b.  In  this  case  the  internal  geometry  is 
a  square  with  a  finite  radius  of  curvature  at  each  corner.  Due  to  the 
close  proximity  of  the  outer  boundary,  an  analytical  solution  is  not 
available.  In  lieu  of  this  a  parametric  study  was  made.  Two 
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parameters  were  necessary  to  describe  the  square:  a/b  the  port  frac¬ 
tion  and  a/p  the  fillet  radius  factor.  The  results  of  this  study 
are  shown  in  Figure  23.  The  limit  points  required  for  this  series  may 
then  be  determined  from  this  curve.  A  typical  photoelastic  picture  of 
this  limit  configuration  is  shown  in  Figure  24. 

The  actual  curves  used  t  >  present  the  data  are  obtained  by  holding 
d/2 p  fixed  and  plotting  the  factor  H  as  a  function  of  a/p  for 
several  values  of  a/b  .  Each  member  in  the  resulting  family  of  curves 
is  located  by  3  to  4  test  points.  In  general  the  correlation  of  these 
tests  in  establishing  a  single  curve  is  excellent.  It  will  be  noted 

that  the  limit  points  are  indicated  by  a  vertical  crosshatched  line 

( 

and  are  a  function  of  d/2p  only.  For  large  values  of  a/b  the  limit 
ing  geometry  cannot  oe  realized,  as  the  limiting  square  is  larger  than 
the  outer  circular  boundary. 

In  general  the  curves  are  quite  smooth  except  near  the  limit  point 
where  a  rapid  increase  in  the  factor  H  is  indicated.  Also  it  will 
be  noted  that  the  factor  H  is  significantly  reduced  for  values  of 
d/2p  >  1  (with  the  exception  of  points  in  the  neighborhood  of  the 
limit,  points).  Also,  as  indicated  in  Figure  25,  a  minimum  value  of  H 
is  obtained  for  fixed  values  of  the  other  remaining  parameters.  The 
value  of  d/2 p  where  this  minimum  occurs  is  not  constant  but  is  a 
function  of  the  other  parameters. 

In  the  simple  slot  tests  for  values  of  N  from  N  ■  3  to  N  ■  8 

-1/3 

the  empirical  rule  H  «s  N  was  observed  (See  Figure  14) .  If  the 
same  rule  were  to  hold  for  d/2 p  j*  1  ,  then  the  results  of  the 


current  series  of  tests  could  be  lnioedlately  extended  for  N  j*  4  . 

To  test  this  hypothesis  a  series  of  six  tests  for  N  *  2,  3,  4,  5,  6, 

and  7  was  conducted  with  the  remaining  three  parameters  held  constant 

at  d/2p  ■  3  ,  a/b  ■  60%  ,  a/p  *  13  .  The  results  of  the  tests 

showing  the  variation  with  N  are  shown  in  Figure  26.  The  dotted  line 

•1  #3 

indicates  the  curve  which  would  be  anticipated  if  an  N  '  rule 
governed  the  behavior.  This  curve  represents  the  data  particularly 
well  in  view  of  the  rather  peculiar  shape  of  the  curve.  The  curve 
shown  has  been  drawn  to  give  the  best  fit  to  all  the  data,  rather  than 
to  tie  to  the  N  ■  4  results.  The  largest  disagreement  of  the  two 
curves  shown  is  6Z. 

It  would  be  advisable  to  make  additional  tests  to  obtain  cross¬ 
plots  similar  to  Figure  26  for  other  fixed  values  of  d/2 p  ,  a/b  , 

-1/3 

and  a/p  to  determine  the  limits  of  applicability  of  an  N  rule. 

This  work  however  was  not  within  the  scope  of  this  investigation. 
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VIII. 


POSITIVE  WEDGE  ANGLE  TESTS 

A  series  of  tests  was  conducted  to  determine  the  effect  of  a 
positive  wedge  angle  on  the  stress  factor  H  .  The  geometry  used  in 
this  series  is  defined  in  Figure  16.  It  is  completely  specified  by 
four  parameters:  N  the  number  of  star  points,  a/b  the  port  frac¬ 
tion,  a/p  the  fillet  radius  factor,  and  QL  the  wadge  angle.  A 
total  of  63  tests  were  conducted  for  N  *  4  ;  a/b  *  40%,  50%,  60%, 

70%,  and  80%  ;  CL  »  10°,  20°,  30°,  and  40°  ;  and  an  appropriate 

range  of  a/p  .  In  addition  five  tests  were  conducted  for  values  of 

N  other  than  4  to  determine  the  effect  of  N 

The  results  for  N  *  4  are  shown  in  Figures  27  through  30.  In 
each  of  these  figures  QL  and  N  have  been  held  constant  and  H  is 
plotted  as  a  function  of  a/p  for  various  values  of  a/b  .  Each 
curve  is  defined  by  three  test  points  and  a  limit  point  at  a/p  *  1 
At  this  limit  point  the  geometry  reduces  to  a  circular  port  for  which 
the  value  of  the  stress  may  be  easily  calculated.  A  typical  isochro- 
matic  fringe  pattern  for  this  geometry  is  shown  in  Figure  31. 

If  these  curves  are  compared  to  the  simple  slot  tests  of  Figure  8, 
which  corresponds  to  this  geometry  with  QL  »  0  ,  it  is  observed 

that  even  moderate  positive  wedge  angles  lead  to  significant  reductions 
in  stress.  This  is  clearly  seen  in  the  crossplots  presented  in  Figures 
32  and  33,  where  H  is  plotted  versus  0L  for  various  values  of  a/b 
with  a/D  held  constant. 
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The  results  of  the  tests  showing  the  variation  of  H  with  N 

are  shown  in  Figure  34.  Again  the  stress  factor  H  decreases  as  N 

-1/3 

increases,  and  the  N  rule  (dotted  curve)  is  in  reasonable  agree¬ 


ment  with  the  data. 


IX. 


NEGATIVE  WEDGE  ANGLE  TESTS 


A  series  of  tests  was  conducted  to  determine  the  effect  of  a  nega¬ 
tive  wedge  angle  on  the  stress  factor  K  .  The  geometry  used  in  this 
series  is  shown  in  Figure  17.  The  definition  of  this  geometry  requires 
four  parameters.  They  are:  N  the  number  of  star  points,  a/b  the 
port  fraction  ,  a/p  the  fillet  radius  factor,  and  the  wedge 
angle.  The  wedge  angle  has  been  denoted  by  ,  a  positive  number,  for 
convenience.  It  will  be  noted  that  this  corresponds  to  an  extension 
of  the  series  reported  in  the  previous  section  for  negative  values  of 

a  . 

Tests  were  conducted  for  N  *  4  ,  ■  401,  5 OX,  6 OX,  707.,  and 

80X  ,  •  10°,  20°,  and  40°  ,  and  an  appropriate  range  of  a/p 

An  additional  series  of  tests  was  conducted  for  values  of  N  from  2 
to  8. 

A  typical  isochromatic  fringe  pattern  is  shown  in  Figure  35.  The 


results  for  N  ■  4  are  shown  in  Figures  36  through  39.  In  each  of 
these  figures,  Q  and  N  have  been  held  constant  and  H  is  plotted 
as  a  function  of  a/p  for  various  values  of  a/b  .  Each  curve  is 
defined  by  three  test  points  and  a  limit  point  at  m/p  ■  1  At  this 
limit  point  the  geometry  reduces  to  a  circular  port  for  which  the  value 
of  the  stress  may  be  easily  calculated.  The  curve  for  0  again 

corresponds  to  a  simple  slot  configuration  (See  Figure  8).  The  curve 
for  ■  30°  has  been  obtained  by  interpolation. 
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The  results  of  the  tests  showing  the  variation  of  H  as  a  func¬ 


tion  of  N  are  «how  i  in  figure  40.  It  will  le  noted  that  the  stress 

.1/3 

decreases  as  K  increases  and  that  the  N  rule  (dotted  lice) 
closely  represents  the  data.  This  is  particularly  true  for  values  of 
N  greater  than  3.  In  the  limit  when  N  ■  2  ,  the  disagreement  is 

only  6.7X.  However  it  must  be  emphasized  that  this  validates  this 
rule  only  within  the  range  tested. 

Cro  splots  of  the  data  are  presented  in  Figures  41  through  43. 

In  these  plots  N  and  a/p  are  held  constant  acd  H  is  plotted  as 
a  function  of  (3  for  various  values  of  a/b  .  The  limit  point  for 


/? 


0  corresponds  to  a  simple  slot  configuration;  the  other 


extreme 


.  a 


m  90°  ,  corresponds  to  a  square  with  finite  radius  of 

curvature  at  the  corners.  Both  of  these  configurations  have  been 
investigated.  It  is  seen  from  these  plats  that  the  variation  of  K 
as  a  function  of  /3  is  relatively  small  for  angles  up  to  50°. 
Particularly  notice  that  for  smaller  port  fractions  the  value  of  H 
Is  essentially  constant  for  /3  <  45°  For  this  reason  the  curve 


for 


13- 


30°  was  obtained  by  interpolation  of  this  data. 
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X. 


ELLIPTICAL  SLOT  TIP  TESTS 

The  nature  of  the  stress  concentration  at  a  star  tip  is  determined 
by  the  star  tip  geometry.  In  the  case  where  the  slot  is  terminated  in 
a  semicircular  tip,  one  point  of  stress  concentration  occurs,  and  it 
is  located  on  the  centerline  of  the  slot  (Figure  2).  On  the  other 
hand,  when  the  slot  is  terminated  with  a  tip  similar  to  those  used  for 
the  slot  width  effect  tests,  two  points  of  stress  concentration  occur 
(Figure  3).  With  these  observations  in  mind,  consider  the  nature  of 
the  stress  concentration  where  the  slot  is  terminated  by  a  semiellipse 
(Figure  18).  If  we  define  6  to  be  the  ratio  of  the  minor  axis  of 
the  ellipse  to  the  major  axis  of  the  ellipse,  then  as  6  — ►  1  the 
star  tip  approaches  a  semicircle  (Figure  44),  and  a  single  point  of 
stress  concentration  would  be  expected.  On  the  other  hand,  when  € 
is  small  (€  — 0),  the  star  tip  resembles  the  type  used  in  the  slot 
width  effect  tests,  and  two  points  of  concentration  would  be  expected. 
It  is  clear  that  at  some  intermediate  value  of  €  the  transition  from 
one  point  of  concentration  to  two  points  of  concentration  must  occur, 
and  for  this  transition  geometry  the  stress  must  be  quite  uniformly 
distributed  around  the  tip.  Intuitively  one  feels  that  this  tuore 
uniform  stress  distribution,  relatively  free  from  concentration  points, 
will  be  a  minimum  stress  configuration. 

To  test  this  hypothesis,  a  series  of  tests  was  conducted  with 
elliptical  star  tip  geometries,  with  values  of  €  ”1,  0.8,  0.6,  0.4, 
and  0.35.  In  all  cases  it  was  found  that  as  €  is  varied  while  all 
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other  parameters  of  the  geometry  are  held  constant,  a  configuration  of 
minimum  stress  occurs  for  a  value  of  £  in  the  range  0.35  <  €  < 
0.70  .  As  anticipated,  this  minimum  stress  configuration  is  found 

at  the  value  of  £  where  the  transition  from  one  concentration  point 
to  two  concentration  points  is  occurring. 

The  geometry  of  a  star  grain  with  elliptical  star  tips  can  be 
specified  by  four  dimensionless  parameters,  a/b  ,  N  ,  a/p  1  , 

and  €  .  The  majority  of  the  tests  were  conducted  for  N  *  4  . 

Typical  fringe  patterns  are  shown  in  Figure  45.  Curves  of  H  (stress 
factor)  versus  a/p  '  ,  for  variou_  /alues  of  a/b  and  €  ,  are 

shown  in  Figures  46,  47,  48,  and  49.  Crossplots  of  H  versus  € 
for  fixed  a/b  and  a/p  '  are  shown  in  Figures  50  and  51.  The 
shaded  area  Indicates  the  value  of  C  at  which  a  minimum  H  occurs. 

In  addition  to  the  tests  for  N  ■  4  ,  several  tests  were  con¬ 
ducted  for  other  values  of  N  .  The  relationship  between  H  and  N 
for  fixed  a / p '  ,  a/b  ,  and  €  is  shown  in  Figures  52  and  53. 

The  dotted  line  indicates  the  variation  of  H  if  the  empirical  rule 
-1/3 

H  ~  N  were  valid.  The  close  agreement  between  the  dotted 

curve  and  the  experimental  values  indicates  that  once  again,  within 

-1/3 

the  range  tested,  the  N  '  rule  adequately  expresses  the  dependence 
jt  H  on  N  .  This  indicates  that  the  extensive  results  obtained 
for  N  *  4  can  be  directly  extended  to  N  +  4  . 

A  series  of  fringe  pattern  pictures  is  shown  in  Figure  45  to 
illustrate  how  the  transition  from  one  point  of  stress  concentration 
to  two  points  of  concentration  occurs  as  €  is  varied.  The  pictures 
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were  taken  with  a  load  of  350  psi,  and  in  each  configuration  N  *  4  , 
a/p  ’  “6.5  ,  and  a/b  *  607«  .  Note  in  particular  the  uniform 
3tress  along  the  star  tip  in  Figure  45-c  where  £  is  close  to  the 
transition  value  as  compared  to  Figure  45-a  or  45-d. 
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XI. 


APPLICATIONS 


Internal  Pressure 


Maximum  Stress  Calculations.  In  the  introductory  section  it  is 

shown  that  the  test  result  can  be  applied  to  the  case  of  internal 

pressure  p|  and  external  pressure  p^  by  superimposing  a  hydrostatic 

stress  state  CT  *  O"  *  O!  =  Ol  *  -p!  .  Thus,  for  this  case, 

xx  yy  1  2  ’ 

•  •  « 

P, 


Po  ’  p< 


°T  =  H(PI  '  Po>  ■  PI 


(17) 


The  factor  H  is  obtained  from  the  parametric  curves  and  is  equal  to 

-<Vpo  • 

If  the  grain  is  case-bonded,  then  to  a  good  degree  of  approxima¬ 
tion^  the  case  can  be  replaced  by  an  equivalent  uniform  external 
pressure  p'  where 


P' 


(18) 


In  this  equation  t£  is  the  thickness  of  the  case  and  and  1/ 

describe  the  elastic  properties  of  the  case.  The  equation  has  been 
derived  for  a  thlckwalled  cylinder  grain.  The  validity  of  this 
equation  becomes  questionable  for  a  star  grain  when  the  web  becomes 
thin.  It  is  clear  that  if  the  web  is  thin  the  external  pressure 
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exerted  by  the  case  will  not  be  uniform,  and  shear  transfer  between 
case  and  grain  may  become  appreciable. 

When  this  equation  is  used  to  find  p'  for  a  star  grain,  a  value 
for  a/b  must  be  chosen  which  properly  represents  the  star  grain  as 
an  "equivalent  thickwalled  cylinder"  (equivalent  in  the  sense  that  p' 
calculated  using  the  chosen  value  of  a/b  is  equal  to  the  average 
value  of  normal  stress  which  the  case  exerts  on  the  actual  star  grain) 
The  determination  of  an  "equivalent  thickwalled  cylinder"  is  a  key 
step  in  the  application  of  photoelastic  test  data  to  case-bonded  star 
grains.  Engineering  intuition  dictates  that  the  equivalent  cylinder 
lies  between  an  equal  web  fraction  cylinder  and  an  equal  area  cylinder 
The  problem  of  choosing  an  equivalent  cylinder  will  be  pursued  in 
detail  in  a  later  section. 

As  a  practical  note  for  calculation  purposes,  it  will  be  found 
that  when  typical  numerical  values  are  substituted  into  equation  (18) 
the  denominator  will  be  approximately  equal  to  1,  indicating  that 
p'  9s  Pj  •  Since  the  stress  calculation  (equation  17)  depends  upon 
the  factor  (p^  -  p'),  p'  must  be  accurately  determined.  Computa¬ 
tional  labor  is  minimized  if  the  denominator  is  put  into  the  form 
1  +  A  .  Then, 


i  -  A  +  A2  +  o(A3)  V 


U9) 


and  an  adequate  expression  for  p'  will  be  obtained  in  most  practical 

cases  by  retaining  the  first  three  terms  in  the  expansion;  i.e., 

/ 
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A  «  A  .  The  value  of  A  is  easily  calculated  from  the 
expression 


A 


Usi _ 

2(1  -  V) 


1  -  2Z/  + 


1  +  V 


1  +  V 


Eb 


Me 


(20) 


In  those  problems  where  V  ~  0.5  and  Ectc  »  Eb  ,  this  approach 
will  be  found  to  significantly  reduce  the  computational  labor.  The 
maximum  stress  (equation  17)  becomes 


^T 


H(  A-  A2)  -  1  Pi  +  0(HPiA3)  (21) 


Maximum  Strain  Calculations.  The  calculation  of  maximum  strain 
follows  from  equations  (2).  For  the  particular  case  of  external 
pressure  p^  and  internal  pressure  p|  ,  which  includes  case-bonded 
grains  as  Indicated  in  the  previous  section,  the  maximum  strain  occurs 
at  the  point  of  maximum  stress  and  is  given  by 


€T  -  k±}L  J  (1  -  I/)H  +  2V  -  lj  p[  +  (  V  -  1)HP;  l  (22) 


In  terms  of  A  for  a  case -bonded  grain,  this  equation  becomes 


—  -  Lv~-  <(1  -  vu  A-  A2)h  -  (1  -  21^) 

Pi  8 


,+  o(BA3) 


(23) 


If  i'1  **  0.5  ,  then  depends  strongly  upon  H  ,  A  ,  and 

]/  .  The  dependence  upon  ]/  is  more  clearly  expressed  by  letting 
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^  -  0.5  - 


4 


S  ,  S  «  i 


Then  equation  (23)  becomes 


!i « i  +  v 

Pi  E 


%<  A  -  A2)  + 


*\ 


It  is  clear  from  this  expression  that  under  these  conditions 
(  J/  sc  0.5),  accurate  strain  calculations  require  very  accurate  values 
of  V  and  interface  pressure  p1  .  It  is  also  clear  that  errors  in 
interface  pressure  are  greatly  magnified  when  H  is  large. 

The  extension  of  these  calculations  to  the  case  of  a  linear 
viscoelastic  material  is  straightforward  as  long  as  the  boundary  con¬ 
ditions  are  of  the  proportional  loading  type.  The  elastic  constants 
E  and  1/  are  replaced  by  appropriate  differential  operators  and 
integrations  carried  out  where  required.  In  particular,  it  should  be 
noted  that  equation  (18)  becomes  a  differential  equation  for  p'(t)  , 

and  the  time  dependence  of  p'(t)  will  in  general  not  be  the  same  as 
pj^(t)  .  Thus,  for  a  viscoelastic  grain  in  an  elastic  case,  the 

assumption  of  proportional  loading  may  not  be  valid.  Treating  the 

> 

problem  as  one  with  proportional  loading  introduces  further  approxi¬ 
mation  into  the  solution. 


Itetwd  adaism 

If  a  case-bonded  solid  propellant  rocket  grain  ia  strain  free  at 
some  reference  temperature,  say  T  •  0  ,  and  if  its  temperature  is 

changed,  strain  will  be  induced  due  to  two  causes.  First,  strain  will 
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be  induced  due  to  the  thermal  expansion  of  the  body.  If  the  body  is 
unrestrained  and  the  temperature  is  uniform  from  point  to  point  in 
the  body,  then  the  body  will  remain  stress  free  although  a  nonzero 
strain  field  exists.  If  however  the  body  is  restrained  in  some  manner, 
these  restraints  will  produce  a  stress  fJeld  which  then  gives  rise  to 
additional  strains  through  the  normal  Hookean  constitutive  equations. 
This  is  expressed  in  the  following  equations. 


«rr  "  i  +  O’,,)]  +  0.^1 

€  -  z  cr  -  o;  +  aijl  +citAt 

ee  e  L  e©  zz  rr  J  T 

€.«*i  o~QQ  )]  +aTAT 


It  , 

G  r© 


It 

G 


©z* 


It 

G  rz 


(25) 


In  each  of  these  expressions  for  normal  strain,  we  see  the  strain 
arises  from  two  sources.  At  the  temperature  change  and  the  Induced 
stresses.  In  these  equations  At  Is  the  temperature  change  and  CL^ 
is  the  coefficient  of  thermal  expansion.  Note  that  the  thermal  portion 
of  the  stress  field  affects  only  the  normal  stress-strain  relations, 
as  thermal  expansion  causes  no  shearing  strains. 

Nov  consider  the  thermal  stress  problem  in  a  long  circular  cylin¬ 
drical  rocket  grain  with  constant  temperature  change  At  ,  a  thin 
elastic  case,  and  a  concentric  circular  internal  port.  The  assumption 
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is  made  Chat  Che  ends  of  grain  are  restrained  in  such  a  manner  Chat 


the  problem  reduces  to  one  of  plane  strain;  i.e.,  €  ■  0  .  The 

equations  for  the  stresses  and  strain*  in  i.he  propellant  are: 


(T 

rr 


^©e 


CT 
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(1  -  a2/b2) 


1  -  a2/r2 


2  2 

(1  -  a  /b 


] 

-|\  f  ,2/r2  j 
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us*  i  r 

-  a2/b2)  L 


(1  -  a  /b  ) 


(1  -  2i</)  -  «2/r2  +  (1  +  1/)CItAt 


(26) 


*  l/,)p'  f(l  -  21/)  +  a2/r2  1  +  (1  +  I/XX-At 

-  *2/b2)  L  J  1 


}  •  (b)2l  [  TTu)  [(1  *  vvxt  -  (1  +  ^c^,]  At 


b)*— 


It  should  be  noted  that  in  the  above  set  of  equations  the  stresses 
and  ^Q0  result  from  the  term  involving  the  interface  pressure 
p'  only.  That  is,  the  tvo-dimenslonal  stress  field  is  completely 
equivalent  to  a  pressure  problem  vith  the  value  of  the  external  pressure 
being  given  by  p’  of  the  above  set  of  equations. 

It  should  be  noted  that  the  above  expression  for  p'  does  not 
agree  with  the  expressions  given  in  Reference  6  and  Reference  9. 
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This  p'  has  been  based  on  the  assumption  6  *  0  where  6  Is 

zz  zz 

the  total  z  strain;  i.e.,  mechanical  strain  plus  thermal  strain. 

The  referenced  reports  give  a  value  for  p'  bused  on  holding  only 
the  mechanical  part  of  6  zero. 

ZZ 

When  the  internal  boundary  is  not  a  circular  port,  it  is  not  possi 
ble  to  write  expressions  for  the  stress  field  similar  to  the  pressure 
problem.  The  assumption  is  made  that  the  effect  of  the  case  may  be 
replaced  by  a  uniform  pressure  given  by  p'  of  equation  (26).  Here 
again  an  equivalent  circular  port  grain  must  be  selected  to  calculate 
p'  .  In  this  case  the  equal  web  fraction  would  yield  a  lower  bound  on 
the  stress  while  ^he  equal  area  equivalence  would  yield  an  upper  bound. 

Strain  calculations  are  most  readily  made  by  first  calculating  v 
the  complete  stress  state  using  equation  (17)  to  calculate  the  maxiqum 
tangential  stress,  the  Internal  bo'tndary  condition  to  calculate  the 
normal  stress,  and  the  condition  of  plane  strain  for  the  axial  stress 
j^i.e.,  €zz  •  0  ,  which  implies  &zz  ■  l/(  +  <7^)  -  eOL^AtJ  . 

Once  the  stresses  are  known,  the  strains  may  be  determined  directly 
from  equation  (25).  It  should  be  noted  that  the  expression  given  for 
maximum  strain  in  equation  (22)  is  incomplete  when  thermal  strains 
are  considered. 

Recently  the  results  of  a  limited  number  of  computer  solutions 
of  case-bonded  and  free  star  grains,  loaded  by  pressure  and  thermal 
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Comparing  these 


(7) 

loads,  have  been  published  by  Becker  and  Brisbane, 
results  with  the  ones  obtained  using  photoelastic  data  and  the  tech¬ 
niques  outlined  in  this  report,  one  finds  remarkable  agreement  for  the 
free  grain  for  both  stress  and  strain.  For  case-bonded  grains,  the 
maximum  stress  calculations  show  good  agreement  when  p'  is  based  on 
an  equal  web  fraction  equivalent.  However,  as  indicated  in  a  previous 
section,  considerable  error  can  occur  in  the  strain  calculations, 
especially  for  large  values  of  H  and  1/  ss  0.5  since  the  strain 
calculation  is  particularly  sensitive  to  errors  in  p1  .  In  the 
examples  of  Becker  and  Brisbane,  V  *  0.4987  and  11  <  H  <  18.2  . 

In  more  typical  configurations  4  <  H  <  10,  so  that  these  examples 
probably  indicate  an  upper  bound  to  the  expected  error. 

Using  the  examples  of  Becker  and  Brisbane  and  the  example  from 
the  Numerical  Example  Section  of  this  report,  tentative  rules  for 
finding  a  better  "equivalent  cylinder"  can  be  derived.  The  examples 
show  that  the  correct  equivalent  is  much  closer  to  equal  web  fraction 
than  to  equal  area,  indicating  that  the  portion  of  the  star  grain 
material  inside  an  equal  web  fraction  cylinder  is  not  effectively 
utilised  in  providing  hoop  stiffness.  The  area  of  this  Ineffective 
material  is  proportional  to  £(a/b)y  -  (a/b)j^  J  where  (a/b)w 
refers  to  an  equal  web  fraction  equivalent  and  (a/b)A  refers  to  an 
equal  area  equivalent.  By  considering  this  area  to  be  201  effective 
in  Increasing  hoop  stiffness,  considerable  improvement  in  the  accuracy 
of  strain  calculation's  is  obta  ined.  In  other  words,  for  purposes  of 
calculating  the  interface  pressure  p*  ,  an  equivalent  cylinder 
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(a/b)_  is  calculated  as  follows: 
E 


(a/b). 


\l 


(a/b -  0.20  «a/b)*  -  (a/b)H 


1 


(27) 


The  stresses  and  strains  calculated  using  (a/b)y  ,  (a/b)^  ,  and 

(a/b)£  are  tabulated  in  Table  II  along  with  the  corresp'-  ding  computer 

solutions.  In  these  examples,  when  (a/b)  .  L\  used  ro  calculate  p' 

the  maximum  error  in  stress  is  less  than  10%  tad  the  naximum  error  in 

strain  is  less  than  15%.  Additional  coroarison  problems  would  be 

useful  to  verify  the  validity  of  the  method.  However,  based  on  the 

information  at  hand,  it  would  appear  that  the  use  of  (r/b)  offers 

E 

considerable  improvement  over  (*/h)„  or  (a/b). 

m  A 


Numerical  Example 

Internal  Pressure.  A  numerical  example  will  be  presented  to 
illustrate  the  dependency  of  the  maximum  t trass  (or  strain)  on  tha 
choice  of  the  "equivalent1*  circular  port  grain.  Assume  that  the  cross 
section  of  the  grain  if  of  the  simple  slot  family  shown  in  Figure  7-a 
with  tha  following  values  of  parameters: 

N  •  4  d/2p  »  1 

a/b  •  66.7%  .a/p  -  8 

Bence  from  Figure  8  the  value  of  the  stress  factor  is  datarmined  to 
be  H  -  9.6  . 
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n.4 


Assume  the  following  properties  for  the  propellant  and  the  case:  J 

f 

%  -  30  x  106  psi  £  «  103 

Uc  -  0.30  l'-  0.49 

tc  “  0.06"  b  -  3"  j 

i 

! 

This  problem  is  identical  to  the  problem  used  to  compare  computer  > 

solutions—obtained  by  ISC  and  by  Rohm  ar.d  Haas.^  Thub  the  maximum  • 

tangential  stress  and  strain  are  known  from  an  independent  source  and  I 

j 

can  be  compared  to  the  results  obtained  using  the  much  simpler  methods  * 

of  this  report .  The  computer  solutions  gJv  .  | 

\ 

3 

v><y/p£  ®  —  0.698  1 

$ 

* 

l 

€  -j/Pi  2.0  x  10”^/psi  j 

In  other  words,  for  an  internal  pressrre  of  100  psi,  the  maximum  tan¬ 
gential  stress  is  *  -  69.8  psi  and  the  maximum  tangential  strain 

is  €rf  -  27,  , 

To  illustrate  the  dependence  of  the  solution  on  the  proper  choice 
of  an  "equivalent  cylinder,"  the  solution  will  be  obtained  using  the 
equal  web  fraction  equivalent,  the  equal  are..  equivalent,  and  the 
effective  area  equivalent  defined  by  equation  (27). 

For  the  equal  web  fraction  equivalent,  (a/b)w  ■  0.667  .  Sub¬ 

stituting  this  value  into  equation  (20),  the  pressure  factor  A  can 
be  calculated  as  follows 


A  -  0.0260 
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In  this  example,  since  A  <<  A  ,  the  use  of  only  the  first  two 
terms  in  the  series  is  completely  justified. 

The  average  interface  pressure  exerted  by  the  case  on  the  grain 
can  be  calculated  using  equation  (19) 

p'  -  Pi(l  -  £a  -  A2^|) 

-  0.9747  pt 

Thus,  for  this  example,  the  internal  pressure  is  transmitted  with  very 
little  loss  to  the  case.  This  is  a  typical  result  for  V  ^  0.5  end 
a  relatively  stiff  case  (i.e.,  Eb/Ectc  <<  1). 

The  maximum  tangential  stress  (tangent  to  the  star  boundary)  is 
calculated  from  equation  (21) 

CTT  -  -  0.757  pj^ 

The  maximum  tangential  strain  is  calculated  from  equation  (24) 

CJj  ■  1.54  x  10”4  p^ 

Thus  the  tangential  stress  is  negative  (compressive)  while  the  tan¬ 
gential  strain  is  positive.  It  should  be  emphasized  that  this  stress 
is  the  -axtmum  (absolute  value)  value  of  tangential  stress  which  occurs 
at  the  star  port  surface.  The  term  tangential  is  used  to  indicate  the 
direction  parallel  to  the  star  boundary  at  the  point  where  the  maximum 
occurs  and  does  not  in  general  mean  the  "0"  direction,  although  in 
this  particular  example  (simple  slot)  the  two  directions  do  coincide. 
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This  maximum  tangential  stress  is  not  necessarily  the  maximum  stress 
occurring  at  the  point..  In  order  to  check  this,  the  other  two  princi¬ 
pal  stresses  must  be  calculated  using  the  boundary  condition  0~n  * 

-  p^  and  equation  (3) 

^  -  >y(crxx+  cryy)  crT) 

■  -  0.86  p^ 

The  three-dimensional  state  of  stress  for  p^  ■  100  psi  is  shown 
‘  elow 

azz-  +  <rT)  -  -  86 


In  this  example  CTT  is  the  smallest  (in  absolute  value)  of  the 
principal  stresses. 

The  same  example  is  now  worked  out  using  the  equal  area  equiva¬ 
lent.  The  port  area  of  the  grain  is  found  to  be 

A  -  3.64  in2 


Thus  an  equal  port  area  circular  port  grain  would  have  a  port  radius 


a  -  1.08 
(a/b)A  -  0.36 

The  pressure  factor  A  is  calculated  from  equation  (20) 

A  -  0.139 

Thus  the  interface  pressure  p'  (equation  19)  is 

p'  ■  0.88  p^ 

Equation  (21)  now  gives 

CFj  -  0.122  Pj^ 

The  strain  is  again  calculated  from  equation  (24) 

€  T  »  3.28  x  10  ^  pi 

Note  that  the  calculation  based  on  equal  port  area  gives  a  very  poor 
estimate  of  both  stress  and  strain. 

As  a  final  example,  the  calculation  will  be  repeated  using  the 
circular  port  equivalent  defined  by  equation  (27) 

(a/b)E  -  0.618 
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The  pressure  factor  A  is  again  calculated  from  equation  (20),  using 
a/b  -  0.618  . 


A  -  0.0336 


From  equation  (21) 

CTT  -  -  0.688  ?i 

and  from  equation  (24) 

€t  -  2.08  x  10'4  pi 

In  this  particular  example,  when  (a/b)^.  is  used  to  calculate  p'  , 
the  error  in  stress  is  1.4%  and  the  error  in  strain  is  4%,  assuming 
that  the  computer  solution  is  correct.  In  all  of  the  comparison 
problems  known  to  the  authors  at  the  present  time,  (a/b)E  yields  the 
best  estimates  of  stress  and  strain  for  a  case -bonded  grain  when  com¬ 
pared  to  the  computer  solution  (Table  II). 


Thermal  Shrinkage.  The  result  of  a  computer  solution  for  the 
thermal  shrinkage  problem  can  be  found  in  Reference  7.  The  configura 
tion  is  a  simple  slot  four  star  grain  with  the  following  parameters : 


a  •  8  in 
b  ■  10  in 
C ■  0.5  in 

3 

E  ■  3  x  10  psi 

l/"  0.4987 

H  ■  18.2  (from 
Figure  8) 


t  -  0.15  in 
c 

Efc  ■  30  x  10^  psi 

K  "  °-3 

aT  -  6.5  x  io'5/°f 
ac  •  6.o  x  io'6/°f 
At  -  -  ioo°f 
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The  computer  solution  gives  the  following  values  for  maximum  tangential 
stress  and  strain 

<JT  -  312  psi 
€t  -  8.15% 

For  comparison  purposes,  this  problem  will  be  worked  out  using 
the  methods  of  this  report.  To  find  (a/b)g  »  note  that  (a/b)^  ■ 

0.8  and  (a/b)^  *  0*318  .  Thup  using  equation  (27) 

(a/b)E  -  0.729 

Using  this  value  in  equation  (26),  the  value  of  p'  is: 

p'  -  -  15.8  psi 

The  maximum  tangential  stress  is  now  calculated  using  equation  (17) 
with  external  pressure  p*  and  internal  pressure  p^  ■  0 

<Tt  ■  285  psi 

This  value  is  in  good  agreement  with  the  computer  solution.  On  the 
other  hand,  when  ( a/b )^  or  (a/b)A  used  th®  calculation,  the 
corresponding  values  of  (7"T  are 

<JT  *  18:’  psi  (equal  web  fraction) 

Cj  *  2,760  pal  (equal  area) 


To  continue  with  the  solution,  CT  and  C  are  now  calculated. 

n  zz 

Since  the  internal  pressure  is  assumed  to  be  zero,  CT  »  0  and 

°zz  -  V<-  +  CTn>  -  EaiAl 

*  163  psi 


The  tangential  strain  can  now  be  calculated  from  equations 

«T-  ^>]+aiAl 

-  7.46  x  10"2 


-  7.467. 


Thus  the  strain,  calculated  using  (a/b)  ,  is  in  error  8.57.  for  this 

£ 

problem,  assuming  that  the  computer  solution  is  the  correct  solution. 
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EMPIRICAL  REPRESENTATION  CF  EXPERIMENTAL  DATA 


As  was  indicated  earlier  in  the  section  on  test  results,  the 
data  obtained  for  the  simple  slot  configuration  of  Figure  7-a  may  be 
represented  by  the  empirical  equation 


h  ■  H'i/yHy  [* + 2  v^p  ]  <28) 

It  would  be  desirable  to  determine  similar  empirical  relation¬ 
ships  to  represent  the  data  for  each  family  tested.  It  has  been  shown 

in  the  previous  sections  that  the  variation  with  N  is  adequately 
-1/3 

described  by  the  N  rule. 

Inspection  of  equation  (28)  indicates  that  if  N  and  (a/b)  are 

held  constant  that  a  plot  of  H  versus  would  yield  straight 

lines.  This  may  be  seen  in  Figure  54. 

For  the  slot  width  effect  tests  it  was  shown  that  the  dependence 
-1/3 

on  N  held  (see  Figure  26);  however,  no  additional  simple  func¬ 
tional  dependency  has  been  observed. 

For  the  positive  wedge  angle  tests  the  N"^^  rule  was  shown 
to  hold  (Figure  34),  and  the  dependence  on  the  fillet  radius  factor 
may  be  seen  in  Figure  55.  It  is  seen  here  that  the  functional  fora 
for  the  empirical  relationship  is : 


H 


-  N-1/3 


cx(a. 


a/b)  + 


c2(a, 


a/b) 


h(a,  */b)  (29) 


indicating  that  the  dependence  on  Va^jO  is  linear*  This  is  als0 
true  for  the  negative  wedge  angle  family  and  the  elliptical  slot  tip 

family  as  may  be  seen  in  Figures  56  and  57. 

Judging  from  the  results  shown  in  the  figures,  it  is  seen  that 

the  functional  dependency  of  H  on  the  parameters  N  and  a  Ip  for 
the  simple  slot,  the  positive  wedge  angle,  the  negative  wedge  angle, 
and  the  elliptical  slot  tip  tests  is  adequately  represented  as  des¬ 
cribed  above.  The  task  then  remains  to  determine  the  unknown  func¬ 
tional  form  cf  equation  (29)  for  the  latter  three  familes.  This  work 
is  presently  underway  and  will  hopefully  be  reported  in  a  later  report 
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XIII. 


CONCLUSIONS 


Parametric  photoelastic  investigations  of  four  families  of 
internally  perforated  rocket  grains  have  been  conducted.  The  results 
are  presented  in  graphical  form  showing  the  stress  factor  H  as  a 
function  of  the  four  geometrical  parameters  necessary  to  describe 
each  family.  Good  correlation  of  these  data  with  previous  work  of 
the  authors  and  analytical  and  numerical  solutions  is  shown. 

Methods  of  applying  the  experimentally  obtained  data  to  other  prob¬ 
lems  of  interest  have  been  discussed.  In  particular  the  problem  of  a 
case -bonded  grain  subjected  to  both  pressure  and  thermal  loads  has  been 
considered.  An  engineering  "rule  of  thumb"  is  suggested  which  allows 
quite  accurate  calculations  of  the  stress  and  strain  problems  considered. 

A  limited  number  of  numerical  solutions  for  problems  similar  to 
those  discussed  here  have  recently  been  published  in  Reference  7.  A 
comparison  of  these  and  all  other  known  numerical  solutions  was  made 
with  those  obtained  using  data  obtained  from  this  study.  In  general 
good  agreement  was  obtained.  In  all  cases  the  disagreement  was  less 
than  10X  in  stress  and  151  in  strain. 

Empirical  formulas  are  derived  in  some  cases  to  represent  the 

variation  of  the  stress  factor  as  a  function  of  the  various  parameters. 

In  all  cases  the  variation  with  the  number  of  star  points  has  b^en 
-1/3 

shown  to  be  N  .  In  all  cases  other  than  the  slot  width  tests 

the  functional  dependency  oa  the  fillet  radius  factor  has  been  shown 
to  be  linear  with  a/|0 
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FIGURE  2 

ISOCHROMATIC  FRINGE  PATTERN 
FOR  TYPICAL  SIMPLE  SLOT  CONFIGURATION 
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FIGURE  3 


ISOCHROMATIC  FRINGE  PATTERN 
FOR  TYPICAL  SLOT  WIDTH  TEST  CONFIGURATION 

,  d/2p  -  2,  a/b  -  70%,  a/p  -  7.98,  H  -  8. 
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FIGURE  4 -a 

EXTERNAL  PRESSURE  VERSUS  FRINGE  ORDER 
Test  #29,  Specimen  AFP-SW-24 
N  -  4,  <3/2 p  -  3,  a/b  -  607.,  a/p  -  6.87 
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FIGURE  4 -b 


EXTERNAL  PRESSURE  VERSUS  FRINGE  ORDER 
Test  #29,  Specimen  AFP-SW-24 
N  -  4,  d/2p  -  3,  a/b  -  607.,  fi/p  -  6.87 
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FIGURE  7 

GEO*’: TRY  OF  LIMIT  CASES 
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SIMPLE  SLOT  TEST 
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FIGURE  9 


EXPERIMENTAL 


SLOTTED  GRAIN 


FIGURE  II 


SLOTTED  GRAIN 


SLOTTED  GRAIN 


:gure  13 


N  ”v  Number  of  Slots 
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Detail  A 


FIGURE  18 

GEOtCTRY  OF  TYPICAL  CROSS  SECTION 
FOR  ELLIPTICAL  SLOT  TIP  TESTS 
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FIGURE  19 


SLOT  WIDTH  EFFECT  TESTS 
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SLOT  WIDTH  EFFECT  TESTS 
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SLOT  WIDTH  EFFECT  TESTS 


FIGURE  23 


SLOT  WIDTH  EFFECT 
LIMIT  POINT  TESTS 


FIGURE  24 

ISOCHROMATIC  FRINGE  PATTERN  FOR  LIMITING  SQUARE 
a/b  -  60%,  a/p  -  8.72,  H  -  19.6 
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FIGURE  25 

EFFECT  OF  SLOT  WIDTH  FACTOR 
N  «  4,  a/b  -  7 OX,  a/p  -  12 
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Number  of  Star  Points 
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FIGURE  31 

ISOCHROMATIC  FRINGE  PATTERN 
OF  TYPICAL  POSITIVE  WEDGE  ANGLE  CONFIGURATION 

N  -  4,  »«  20°,  s/b  -  7 OX,  «/p  »  15.12 
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FIGURE  32 


POSITIVE  WEDGE  ANGLE  TESTS 


N  -  4 


a/p  -  8 


FIGURE  33 

POSITIVE  WEDGE  ANGLE  TESTS 


Number  of  Star  Points 


FIGURE  35 


N  * 


(NOTE: 


ISOCHROMATIC  FRINGE  PA1TERW 
OF  TYPICAL  NEGATIVE  WEDGE  ANGLE  TONFIGUP 'TION 

,  p  -  20°,  p/b  -  58% ,  a/p  -  6.33,  H  «■  7.22,  pq  -  300  psi 


Distortion  of  fringe  pattern  on  bottom  extremity  of  model 
was  due  to  de I.iminat ion  of  lens  in  optical  bench.) 


NEGATIVE  HEDGE  ANGLE  TESTS 


NEGATIVE  WEDGE  ANGLE  TESTS 


Experimental 
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«  -  0.6,  2 p  -  0.753 


«  “  0.4,  2p  -  0.753 


FIGURE  44 

DETAIL  OF  ELLIPTICAL  SLOT  TIP  GEOMETRY 
FOR  VARIOUS  ECCENTRICITIES 
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FIGURE  45 

VARIATION  OF  ISOCHROMATIC  FRINGE  PATTERN  WITH  ECCENTRICITY 
pQ  -  350  p*i,  «/b  -  601,  2p  »  0.5" 
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ELLIPTICAL  SLOT  TIP  TESTS 


ELLIPTICAL  SLOT  TIP  TESTS 


FIGURE  48 


FIGURE  49 


8’0 


FIGURE  51 


ELLIPTICAL  SLOT  TIP  TEST 
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In  the  present  investigation  the  study  has  been  extended  to  include  families 
requiring  four  parameters  to  complete  their  description.  Four  additional 
families  have  been  investigated  to  determine  the  effect  of  a)  slot  width,  b) 
positive  wedge  angles  of  tha  star  slot,  c)  negative  wedge  angles  of  the  star 
slot,  and  d)  eccentricity  of  elllptically  shaped  star  tips.  In  each  of  the 
studies  it  it  shown  that  the  N’l/'  rule  holds  approximately ,  and  empirical 
formulas  similar  to  the  one  given  above  have  been  derived  in  some  cases. 

A  taction  on  the  application  of  the  test  results  to  more  general  engineering 
problem*  is  presented.  In  this  section  a  comparison  la  made  between  a  recommended 
method  of  Application  and  a  limited  number  of  numerical  tolutiors  obtained  via 
computer  analysis.  Undar  the  worst  sat  of  conditions  available,  the  strea* 
calculations  agreed  within  10%  and  the  strain  within  15%. 
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